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Abstract. Protostellar jets most probably originate in turbulent accretion disks surrounding young stellar objects. We investi- 
gate the evolution of a disk wind into a collimated jet under the influence of magnetic diffusivity, assuming that the turbulent 
pattern in the disk will also enter the disk corona and the jet. Using the ZEUS-3D code in the axisymmetry option we solve 
the time-dependent resistive MHD equations for a model setup of a central star surrounded by an accretion disk. The disk is 
taken as a time-independent boundary condition for the mass flow rate and the magnetic flux distribution. We derive analytical 
estimates for the magnitude of magnetic diffusion in a protostellar jet connecting our results to earlier work in the limit of ideal 
MHD. We find that the diffusive jets propagate slower into the ambient medium, most probably due to the lower mass flow rate 
in the axial direction. Close to the star we find that a quasi stationary state evolves after several hundred (weak diffusion) or 
thousand (strong diffusion) disk rotations. Magnetic diffusivity affects the protostellar jet structure as follows. The jet poloidal 
magnetic field becomes de-collimated. The jet velocity increases with increasing diffusivity, while the degree of collimation 
for the hydrodynamic flow remains more or less the same. We suggest that the mass flux is a proper tracer for the degree of jet 
collimation and find indications of a critical value for the magnetic diffusivity above which the jet collimation is only weak. We 
finally develop a self-consistent picture in which all these effects can be explained in the framework of the Lorentz force. 
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1. Introduction 



MHD - ISM: jets and outflows - stars: mass loss - stars: pre-main sequence 



Observations of young stellar objects (YSOs) have revealed 
two main features during the phase of star formation - the 
presence of accretion disks and energetic outflows, often ob- 
served as bipolar jets (Mundt et al. 1984, Lada 1985, Ray et al. 
|1996 ). Images and spectra show that these flows are of high- 
velocity (~ 300kms _1 ) and well collimated (opening angle 
< 10° on scales of 1000 AU). The data also suggest that the 
jet collimation must be achieved already close to the central 
source, at distances <100 AU. The jet mass outflow rates are 
~ 10~ 10 - lO _8 M /yr and typically a factor 10-100 smaller 
than the disk accretion rates. 

Besides the special case of protostellar jets, all astrophysi- 
cal jets detected so far seem to be attached to objects where an 
accretion disk is indicated to be present. In particular, this holds 
for jets observed in radio loud active galactic nuclei and quasars 
(Bridle & Perley 1982), highly energetic galactic objects as 
Sco X-l (Padman et al. 1991 >, and microquasars (Mirabel et 
al.|l992|). 



and momentum fluxes exclude the possibility of a thermally 
or radiation pressure driven wind. The observed fluxes are 



much higher than the protostar could provide (DeCampli 198 1 



Konigl 1980- The conclusion is that it is the magnetic field 



which is responsible for protostellar jet formation - accelera- 
tion and co llima tion of the i nitial s tellar or disk wind (Pudritz 
& Norman 1983 , Camenzind 199C ). This magnetic field can be 
generated by some dynamo process either in the central young 
star itself, or the surrounding accretion disk, or it could be pro- 
vided by the interstellar medium as a "fossil" field. What kind 
of mechanism turns the in-flowing matter of the accretion disk 
into an outflow from the disk (or the star) is still not really 
known, although it se ems to be clear that magnetic fields play 
a major role (Ferreira 1997 ). 

Examples of time-dependent MHD simulations of the jet 
formation include models of collimating disk winds, where this 
disk is taken as a boundary condition for the outflow (Ouyed 
& Pudritz 1997, hereafter OP97), models that consider the in- 



teraction of a stellar dipolar magnetosphere with the accretion 



Therefore, the similarities between jets from the differ- 
ent sources imply that the basic mechanism for jet formation 
should be the same. For protostellar jets, the observed mass 



disk as well as the disk structure itself (Miller & Stone 1997 1, or 



some combination of both approaches (Fendt & Elstner 1999 



2000, hereafter FE00). 
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In this paper, we are interested in a time-dependent simu- 
lation of the resistive MHD equations for the corona above the 
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accretion disk around a young star. The underlying disk pro- 
vides a fixed boundary. The presence of a large scale poloidal 
magnetic field provided by the disk is assumed. After a brief 
theoretical introduction in §2, we describe our model and the 
numerical approach in §3, with emphasis on the effect of mag- 
netic diffusion. For the numerical approach, we introduced 
magnetic diffusion in the ZEUS-3D MHD code. Tests of our 
code are included in Appendix. The results for protostellar jets 
are discussed in §4. We compare our diffusive MHD simula- 
tions of jet formation with the non-diffusive case. 

2. Magnetic jets from accretion disks 

Despite a tremendous amount of work conce rning the forma- 
tion of magne tic jets (e.g. Bl andfo rd & Payne 1982 Pud ritz & 
Norman 1983 , Lovel ace et al. 1987 , H eyvae rts & Norman 
Pelletier & Pudri tz |1992| , Li et al. |1992|, Contopo ulos 
Fendt et al. I1995L Kudoh & Shibata 119974 |l997bL Fendt & 



1989 



Memola 2001 ) the mechanism which actually launches the jet 
from the disk remains unclear. Most of the papers dealing with 
the theory of magnetized accretion disks driving jets have been 
following the principal approach of Blandford & Payne as- 
suming stationarity, axisymmetry and self-similarity (Konigl 
1989| , Wardle & Konigl |1993| , Li |1995[ Ferreira |1997| ). After 



all, it is clear now that the launching of a jet from the accreting 
disk can be described as a purely magnetic process. Ferreira 
(1997) has derived trans-Alfvenic, stationary self-similar jet 
solutions with a smooth transition from a resistive disk, where 
the Lorentz forces are actually responsible for lifting the ac- 
creting gas in vertical direction. 

Progress has been achieved, too, in the simulation of the 
time-dependent MHD jet formation from accretion disks. Two 
major, distinct approaches in order to deal with the complex- 
ity of the jet formation process have been undertaken so far. 
One approach is to take the rotating disk as a fixed bound- 
ary condition for the simulation of the jet. Depending on the 
choice of the initial setup (magnetic field, density and pres- 
sure) and the choice of the gravitational potential, the numeri- 
cal results differ in the degree of collimation and the velocity 
of the resulting jet flow. From an initially split-monopole mag- 
netic field configuration, colli mated, non-stationary outflows 
were obtained (Ustyugova et al. 1995 ). In this simulation non- 
equilibrium initial conditions and a softened gravitational po- 
tential were used. For the same c onfigu ration, but for a stronger 
magnetic field, Romanova et al. ( 1997 ) again obtained station- 
ary but only weakly collimated flows. Ouyed & Pudritz ( 1997 ) 
studied the jet formation embedded in a disk corona initially in 
hydrostatic equilibrium and in pressure balance with the disk 
surface. Essentially, after 400 disk rotations, a stationary col- 
limated jet flow emerges. Similar results were obtained by an- 
other recent study taking into account also the time-dependent 
behavior of the disk boundary condition (namely the field in- 
clination) due to the evolution of the disk wind (Krasnopolsky 
et al jl999| ). The main advantage of the approach of a fixed disk 
boundary condition is the numerical stability of simulation over 
a long time scale. The jet launching itself - the process of di- 
verting accreting matter in the disk into an outflow - cannot be 
treated this way. 



The other approach was therefore to include the simulation 
of the disk structure in the simulation. The first step in this di- 
rection was made by Uchida & Shibata (1985) and Shibata & 
Uchida (1985 , 198£| ) in their pioneering work considering time- 
dependent jet formation. Essentially, the authors show that the 
magnetic twist of the magnetic field induced by the rotation 
of the disk gives rise to Lorentz forces pushing the disk ma- 
terial upwards. The back-reaction of the magnetic field on the 
disk (magnetic braking) may lead to a sub-Keplerian disk ro- 
tation. These results have been confirmed by Stone & Norman 
(1994). MHD simulations considering the diffusive accretion 



disk in interaction with a stellar dipolar magnetosphere reveal 
the c ollapse of the inner di sk aft er a few rotations (Hay ashi 



et al. 1996, Goodson et al. 1997, Miller & Stone 1997). At 



the same time, episodic ejection of plasmoids are generated 
in outer parts of the disk wind. A collimated axial jet feature 
is observed. However, probably because of numerical reasons, 
these simulations could be performed only for a few or tens of 
Keplerian periods of the inner disk. Therefore, the results may 
depend strongly on the initial setup. The assumption of con- 
stant diffusivity in the disk and the corona is probably not very 
realistic. 

The long-term evolution of such models has been inves- 
tigated by FE00, however without treating the disk structure 
in the simulation. They find that the axial jet feature observed 
by Goodson et al. (1997) disappears on longer time-scales. A 
two-component quasi-stationary outflow (from disk and star) 
evolves after thousands of rotational periods. This flow is un- 
collimated on the spatial scales investigated, in agreement with 
the observations indicating a jet radius 100 times larger then 
the grid size of the numerical simulations discussed above. 

Recent MHD simulations of the jet formation from accre- 
tion disks by Matsumoto et al. (119961), Kudoh et al. (11998, 



2002b, 2002a ), and Kato et al. (2002) investigate the disk 



jet interrelation and may probably explain the time-dependent 
ejection mechanism in the jets. However, as we have pointed 
out above for the model setup of a magnetic dipole surrounded 
by a disk, these simulations also were carried out for several 
rotation periods only. The question remains as to how the sys- 
tem under consideration behaves on a long time scale. Further, 
the simulations apply the approach of ideal MHD, an approxi- 
mation, which is most probably not strictly valid for accretion 
disks, especially for protostell ar accr etion disks. In this respect, 
the work by Kuwabara et al. (2000) is of particular interest as 
the authors extend the ideal MHD approach and include re- 
sistivity in the jet formation process. Comparing their simula- 
tions to the close environment of a supermassive black hole, 
Kuwabara et al. derive a critical value for the strength of mag- 
netic diffusion. A normalized magnetic diffusivity below 0.05 
may explain the observed activity in active galactic nuclei. In 
this case, the mass accretion and jet launching takes place inter- 
mittently. The paper also demonstrates the difficulty of carry- 
ing out such simulations, as gravity has been treated applying a 
softened potential. On the other hand, this simplification allows 
for a simulation lasting up to 5 rotations of the accretion torus. 

Again, we note that it is just the limitation in the time evo- 
lution which lead us to the decision to take into account the 
accretion disk only as a boundary condition. In this paper, we 
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are mainly interested in the jet formation process (acceleration 
and collimation) and not in the jet launching mechanism from 
the accretion disk. In particular, we think it is essential to inves- 
tigate whether a jet actually survives also on long time scales. 

3. Model setup and equations 

The paradigm accepted in this paper is the jet launched from 
a diffusive, turbulent accretion disk around a young stellar ob- 
ject. We may expect that the turbulence pattern in the disk may 
also enter the disk corona and the jet, and that the jet flow it- 
self is subject to turbulent diffusion. We will discuss this idea 
of a diffusive protostellar jet below. Despite the fact that we 
take into account the effect of magnetic diffusion for the jet 
formation, our model setup is similar to the models in OP97 
and FEOO. 



3. 1. Resistive MHD equations 

In order to model the time-dependent evolution of jet forma- 
tion, the set of resistive MHD equations to be solved is 



| + V-(H = 



at + 



(i) 



V(p + pa)+pV*-^-^=0 (2) 



dB 

~dt 



V x [v x B 



4tt 



-VJ 
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+ (p + p x )(V • v) - -yVJ = 



V ■ B = 



4tt 
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(4) 



(5) 



(6) 



For the gas law we apply a polytropic equation of state, p = 
Kp 1 with a polytropic index 7 = 5/3. Hence, we do not solve 
the energy equation (4). Instead, the internal energy of the sys- 
tem is defined with e = p/(7 — 1). Such a simplification is 
not expected to affect the result of our simulations much, as 
the resistive dissipation term in the energy equation is negligi- 
ble compared to the other terms because of the factor 1/c 2 (see 
also Miller & Stone 1997 ). The magnetic diffusivity is denoted 
by the variable 77 (see Sect. 3.3). Additional to the hydrostatic 
pressure p, an Alfvenic turbulent pressure p\ = p/Pt with 
/3t = const, (see OP97, FEOO) is included. Alfven waves from 
the highly turbulent accretion disk are expected to propagate 
into the disk corona, providing the perturbations for some de- 
gree of turbulent motion also in the jet. As shown by OP97 the 
additional Alfvenic turbulent pressure is able to support a cold 
corona above a protostellar accretion disk as suggested by the 
observations. 

We solve the equations above using the ZEUS-3D code 



(Stone & Norman 1992 1992) in the 2D-axisymmetry option 



for cylindrical coordinates (r, (j), z). We apply a point mass 
gravitational potential $ = — 1/ V r 2 + z 2 located in the origin. 




Fig. 1. Initial setup for the jet simulation. Shown is the part of 
the computational box close to the origin (the part of the "inner 
jet"). The initial hydrostatic density distribution is indicated by 
the thin concentric isocontours. Thick lines denote the initial 
poloidal field lines of a force-free potential field. 

A finite physical magnetic resistivity is added to the original 
ZEUS-3D ideal MHD code. Tests of our now diffusive ZEUS 
code are presented in the Appendix. 

For our computations we normalize the variables to their 
value measured at the inner disk radius n (see OP97, FEOO), 
e.g. p — ► p/pi- The subscript T assigns that the variables 
are taken at this radius. The time is measured in units of a 
Keplerian rotation at the inner disk radius. The normalized 
equation of motion eventually being solved with the code is 



- + (.'.V'X = 

The coefficients /3, = 
Keplerian speed vk,% 



2 f x B' 



V'$'. (7) 



5% Pi p' Si p' 

8-rrpi/B 2 and Si = PiV 2 K l /p l with the 

= \J GM/ri, correspond to the plasma 
beta and the Mach number of the gas at the inner disk ra- 
dius. For a "cold" corona with p' A > 0, it follows (3t = 
l/(5j(7 — l)/7 — 1). In the following we will omit the primes 
and will discuss only normalized variables if not explicitly de- 
clared otherwise. 



3.2. Initial and boundary conditions 

As the numerical simulation of a magnetized disk is still a diffi- 
cult task and is not yet fully resolved, we chose to study the for- 
mation of the jet flow independent of the evolution of the accre- 
tion disk. A precise disk-wind theory would predict the amount 
of angular momentum and energy carried away from the disk. 
Here, we prescribe the disk as a fixed, time-independent bound- 
ary condition for the jet. 

One may question the combination of a magnetically dif- 
fusive disk corona and a steady-state disk magnetic field dis- 
tribution. Naturally, the time scales of the disk evolution are 
definitely longer than those for the jet flow (except for the un- 
known mechanism responsible for the jet knots). Also, for a 
disk jet, the disk is acting as a source for magnetic flux. This 
can be achieved either by a dynamo process working in the 
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disk generating the magnetic field or just by advection of the 
interstellar magnetic field by the disk towards the central star. 
The time scale for both processes is longer than the jet time 
scale and, thus, we may safely assume a fixed magnetic flux as 
accretion disk boundary condition. 

As initial condition we choose the same setup as OP97. The 
initial poloidal magnetic field is defined by the current-free po- 
tential field configuration of the ^-component of the vector po- 
tential, Atf, = (^\/r 2 + (zd + z) 2 — (zd + z)j jr. The dimen- 

sionless disk thickness z<\ satisfying (zd + z) > for z < is 
introduced in order to avoid any kinks in the field distribution. 
The initial coronal density distribution is in hydrostatic equi- 



librium, p = (r 



-3/4 



The initial corona is defined by 



two free parameters Si and 

The disk itself as a boundary condition for the jet flow is 
in centrifugal balance and penetrated by a force-free magnetic 
field. As the disk boundary condition is time-independent, the 
initial potential field magnetic flux from the disk is concerned. 
The toroidal component of the magnetic field in the ghost zones 
(z < 0) is chosen as B^z < 0) = Pi/r, where pi is another 
free parameter. The mass flow rate from the disk surface into 
the corona is defined by the injection velocity and the density 
of the injected material. With the launching angle Oo(r, z = 0) 
(measured from the jet axis), the velocity field in the ghost zone 
is v = (v r ,Vj,,v z ) = Ui n j(vpSin9o,WK,fpCos6o) for r > 1 
with v m j as a free parameter. For r < 1 the inflow velocity is 
set to zero, which actually defines the inner edge of the disk. 
The inflow density is given as p^ = Kjr -3 / 2 , with K\ as a free 
parameter. 

Besides the disk "inflow" boundary condition, the bound- 
ary condition along the symmetry axis is set as "reflecting", 
and along the tw o remaining boundaries as "outflow" (see also 
Stone & Norman [l99^ , |l 992j OP97; FE00). Figure 1 shows the 
initial setup of a hydrostatic density distribution together with 
the potential magnetic field for the part of the computational 
box close to the origin (the region of the "inner jet"). 

Our choice for the free parameters is the following. We 
have <5i = 100, m = 100, p { = -1.0, and u inj = 0.001, 
similar to OP97 and FE00 in order to allow for a comparison 
of the results. For the plasma-/? we choose /3; = 0.282 (which 
is equivalent to OP97, but is based on the original ZEUS-3D 
normalization) or a lower value fii = 0.141 which has some 
numerical advantages [] The lower /3j does not change the gen- 
eral behavior of the jet. The jet evolution is faster (in physical 
time, not in CPU time) and the Alfven surface is slightly shifted 
in z-direction, but the jet internal structure remains very simi- 
lar. 



3.3. Magnetic diffusion in jets 

Most of the models of MHD jet formation deal with the colli- 
mation and acceleration of a plasma flow in the case of ideal 
MHD. However, it seems to us quite possible that magnetic 
diffusivity may play an important role in protostellar jet forma- 



tion. There are (at least) two reasons which may account for 
that. 

The first reason may be the fact that the jet material of 
young stellar objects is not fully ionized, in difference from the 
case of relativistic jets in AGN or microquasars. The fraction 
of ionization (the ratio of ion to neutral particle density) de- 
rived from optical observations is about 0.5 — 0.01 with the ten- 
dency to decrease along the jet (Hartigan et al.1994, Bacciotti 
& Eisloffel 1999). From this it can be expected that diffusive 
effects may indeed affect the MHD configuration. Theoretical 
studies on this topic have just started recently. We refer to Frank 
et al. (1 999 ) who investigated the asymptotic MHD jet equilib- 
ria under the influence of ambipolar diffusion, showing that the 
initial MHD configuration of the jet changes due to ambipo- 
lar diffus ion at least on the parsec scale. Another reference is 
Ferreira ( 1997 ) who showed the essential role magnetic dif- 
fusion plays for the launching mechanism of the jet from the 
accretion disk. 

The other reason for the existence of turbulence in the jet 
formation region is the fact that the jet launching object itself 
- the accretion disk - is highly turbulent. While turbulence is 
an intrinsic property of accretion disks (and necessary for the 
accretion process itself), turbulence can further be driven in 
the disk corona by the differential rotation of the disk, which 
winds-up the coronal magnetic loops leading to powerful re- 
connection processes (see Miller & Stone 2000). 

It seems to be natural to expect that the turbulent pattern 
being definitely present in the disk-jet interaction region be- 
comes also advected with the jet/wind flow into the domains at 
greater height above the disk. We just note that additionally the 
interaction of the jet with the ambient medium leading to vari- 
ous kinds of instabilities increase the turbulent pattern in the jet 
(however, we do not expect that the latter process does affect 
the jet collimation region). 

In Faraday's law (3), the ratio of the first to the sec- 
ond terms in the brackets is the magnetic Reynolds number. 
It can be represented by R m = vL/r], where L is a typi- 
cal length scale and v is a typical velocity. Due to the large 
length scales in the astrophysical context, R m based on the 
microscopic diffusivity is very large compared to unity. For a 
fully ionized hydrogen plasma the microscopic diffusivity is 

is the classical 



7K 



1 See also Appendix A in FE00 for the different choice of /3i in 
OP97 (equal to unity) and here and in FE00. 



~ r c(v t h/c) , where r e 
electron radius, and vth = y/k-B,T/m c is the electron thermal 
speed. For T = 10 4 K, v = lOOkm/s, and L = 100 AU we 
obtain R m = vLjr\^K. 10 15 . 

It is clear that the relevant diffusivity in the protostellar disk 
and jet is most probably "anomalous" determined by macro- 
scopic MHD instability phenomena (see above), with the re- 
sulting magnetic Reynolds number being much smaller. This 
magnetic turbulence we may parameterize the same way as in 
a Shakura-Sunyaev model for the hydrodynamical viscosity. 
We can define a turbulent magnetic diffusivity t/t = a m vL, 
where a m < 1. As a characteristic velocity for the dynami- 
cal change of the system we might choose the poloidal Alfven 
speed, va = B p 



'Anp. If we choose the size of our com- 
putational box ?' max as a typical length scale and a m « 0.1, 
we obtain R m = VA^max/^T ~ 10 as the typical magnetic 
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Fig. 2. Global evolution of the jet on a grid of (z x r)=(280 x 40)r; with a resolution of 900 x 200 elements. Shown is the state of 
evolution after t = 400 rotations of the disk inner radius for different magnetic diffusivity, rj = (top), r\ = 0.01 (middle), and 
r\ = 0.1 (bottom). Thin lines denote 30 logarithmically spaced isocontours of density. Thick lines denote 20 and twenty linearly 
spaced magnetic flux surfaces (or poloidal field lines). The parameters are Si = 100, k; = 100, /i; = —1.0, /?; = 0.282, and 
v inj = 0.001. Note the preserved initial hydrostatic density and force-free field distribution in front of the bow shocks. The figure 
demonstrates that the bow shock advances slower with increasing diffusivity. 



Reynolds number for the global jet evolution. In difference, 
typical length and time scales are different if we are interested 
in the evolution of the local jet structure (as needed for example 
for the definition of the numerical time stepping in the code). 

It might be expected that the diffusivity throughout the jet 
and the disk corona differ considerably. However, it seems nat- 



ural to expect the diffusivity in a corona close to the disk sur- 
face not to differ much from the value in the outer part of the 
disk. For simplicity, and since we are first interested in the gen- 
eral effect, our simulations are performed with a constant dif- 
fusivity parameter. To introduce a non-constant diffusivity is 
straight forward. 
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Fig. 3. The evolution of the inner jet approaches the quasi- 
stationary state. Shown are poloidal magnetic field lines in the 
case of r\ = 0.1 for different time steps, t = 250, 300, 350, 400 
{thick solid, thin solid, dashed and dotted lines). Same param- 
eter setup as in Fig. 2, except /3j = 0.141. Grid size 280 x 80 
elements for a physical size of (140 x 40)r;. The picture shows 
how the poloidal magnetic field lines diffuse outwards but ap- 
proach a (quasi)-stationary state after 400 rotations (see the 
dashed and dotted lines almost coinciding). 





0.5 



Fig. 4. Magnetic diffusivity and global time scales. Left Global 
magnetic Reynolds number as defined from the scale of the 
computational box of the inner jet L = 20, v = va- Shown 
is the i? m -profile along slices in z-direction at r = 7, 13, 20 
(solid, dotted, dashed line, respectively), for the simulation run 
with i] = 0.1, Pi = 0.141. The plot shows the variation of the 
typical velocity v, for which we have chosen the local Alfven 
speed. Right Time of stationarity. This plot shows for different 
magnetic diffusivity the time period when the inner jet reaches 
a quasi-stationary state. This time is estimated from the evo- 
lution of the poloidal magnetic field lines (see Fig. 3) and the 
error bars indicate our uncertainty. 



The local magnetic Reynolds number is also described by 
the ratio of the dynamic to the diffusive time scale R m = 
Tdiff/Tdyn- Here we may define Tdiff = mm(l 2 /ri) and Td yn = 
min(7/vA), where I is the size of the numerical grid cell. 
Because of the internal structure of the jet, these minimum val- 
ues can sometimes be well below the actual global magnetic 
Reynolds number described above, which is determined by the 
characteristic length scale of our simulation L = r max . 



It is interesting to note that in the model of OP97 a turbu- 
lent Alfven pressure, pa, has been introduced in order to estab- 
lish a pressure equilibrium of a cold corona (or jet) above the 
(hot) accretion disk. Although this turbulence effect was taken 
into account for the pressure balance, OP97 did not consider it 
as a reason for a turbulent magnetic diffusivity. However, it is 
straight forward to relate the Alfvenic turbulent pressure pa to 
a turbulent magnetic diffusivity t]t- In order to derive such a re- 
lation, we now choose the turbulent velocity field vt instead of 
the local Alfven speed as the typical velocity for the turbulent 
diffusivity, 



Vt = a m v T L. 



(8) 



With the definition /3t = (c s /vt) 2 and taking into account the 
fact that c? s = ~fp/p for an adiabatic or poly tropic gas law, we 
obtain 



v 2 = J_P 

Vt Pt P ' 

The normalization gives 



1 



p v 



K,i 



or 



/7-1 _ 



(9) 



(10) 



With the condition of sub-Alfvenic turbulence OP97 derived 
Pt = 0.03. With a typical value for Si ~ 100 and a "mean" 
value for the normalized density p 1 ~ 10~ 2 we obtain for the 
normalized magnetic diffusivity 



7/ ~ 0.015 



0.1 



L0 



(11) 



Since the diffusivity changes only weakly with the density, rj ~ 
p 1 / 3 , this provides a good estimate on the strength of magnetic 
diffusion. A self-consistent simulation would take into account 
the relation between diffusion and density as in Eq. (10). For 
comparison, we run a few of simulations with such an setup. So 
far, we find no significant difference to the computations with 
rj = const. 

4. Results and discussion 

We now discuss the results of our numerical simulations con- 
sidering the MHD jet formation under the influence of mag- 
netic diffusion. For this, we have run simulations with a dif- 
ferent parameter setup. The simulations were performed (i) in 
domains of different physical size in order to investigate the in- 
fluence of boundaries and to obtain information about the large- 
scale flow, and also (ii) with different numerical resolution. 

We detected numerical artifacts (a spurious velocity pat- 
tern) in the corners of the grid where outflow boundary condi- 
tions meets the other (inflow, reflecting) boundary conditions. 
In general, this artifacts remain localized close to the corners 
of the grid over many hundred of disk rotations. All the re- 
sults discussed in this paper are not affected by these effects 
as we mainly concentrate on the inner part of the jet flow. We 
performed one reference set of global simulations with high 
resolution (numerical mesh of 900 x 200 grid points, physical 
grid of (z x r) = (280 x 40)r;). In order to investigate effects 
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which concern only the gross behavior the jet flow and not its 
structure in detail, we run another set of simulations with lower 
resolution (numerical mesh of 280 x 80 grid points, physical 
grid of (z x r) = (140 x i0)r{). The much faster computation 
of the low resolution simulations allowed us to follow the jet 
evolution for a very long time even in the case of a high mag- 
netic diffusivity (up to 4000 disk rotations). The computational 
domain was a factor 2 shorter in direction of propagation. All 
other parameters were the same as for the high resolution runs. 

4. 1. Formation of the global jet 

The first point to be noted when comparing the large-scale evo- 
lution of the diffusive with the non-diffusive flow, is the differ- 
ent speed of the bow shock. With increasing magnetic diffusiv- 
ity, the bow shock propagates slower. The bow shock pattern 
velocity is about 0.38, 0.35, 0.28 for r\ = 0, 0.01, 0.1, respec- 
tively (Fig. 2). We will later see that, in apparent contrast, the 
jet velocity increases with 77. 

As in the case of ideal MHD simulations, before the bow 
shock front builds up, torsional Alfven waves propagate from 
the disk surface into the corona slightly modifying the initial 
hydrostatic equilibrium. The super-Alfvenic flow catches up 
with and passes this wave front. 

As the bow shock propagates through the corona it leaves 
behind a cavity of matter with dilute density and high veloc- 
ity. The initially purely poloidal magnetic field becomes more 
and more helical. The toroidal magnetic field component is first 
generated by winding-up the initial poloidal field due to differ- 
ential rotation between the disk and the static corona but later 
comes out as a natural result of the MHD flow due to the iner- 
tial forces of the matter. 

The internal structure of the jet behind bow shock layer is 
smoother in the case of a non-vanishing magnetic diffusion (see 
Fig. 2). We note that also the "wiggly" structure in the density 
distribution close to the disk in the case of 77 = is less promi- 
nent in the case of diffusive simulations. These "wiggles" seem 
to be a numerical artifact probably due to the density jump be- 
tween the disk and the jet, however, the density variations are 
only of the order of some percent. 

Here, we note another important point. In this paper we 
were interested in the cases of a typical MHD jet flow starting 
as a sub-Alfvenic (but super-slow magnetosonic) flow from the 
disk surface, being accelerated to super Alfvenic and super- 
fast magnetosonic sp eed, a s e.g. described in the paper by 
Blandford & Payne ( 1982 ). This case often described as a 
magneto-centrifugally driven disk wind/jet differs from the 
case where matter is injected into the disk corona already with 
super Alfvenic speed. The latter case applies in the case of a 
relatively weak disk poloidal magnetic field as for example in 
the case of a central dipolar field with a strong gradient in ra- 
dial direction. These winds are initially driven by the (toroidal) 
mag netic fi eld pressure g radien t in vertical dire ction ( Lovelace 
et al. |1987| , Contopoulos |l994| , Fendt & Elstner |2000| ). 

In this paper we are interested only in the classical case of 
MHD jet formation. However, in our simulations we note that 
the Alfven surface moves as a function of time until the quasi- 




Fig. 5. The de-collimation of the stationary state poloidal mag- 
netic field due to magnetic diffusivity. Same parameter setup as 
in Fig. 3. Shown is the poloidal magnetic field line distribution 
of the inner jet in the state of quasi-stationarity (t = 400) for 
vanishing diffusivity 77 = {thick) and for 7/ = 0.1 {thin). 



stationary state (see below) is reached. For a smaller than mod- 
erate magnetic diffusivity (77 = 0, ...,0.5), the location of the 
Alfven surface is always within the active zones well above the 
accretion disk boundary. For higher diffusivity the Alfven sur- 
face may advance into the disk for small radii and the character 
of the MHD flow is changed. We do not consider these cases in 
our discussion. 

4.2. Inner jet 

In the following we discuss the evolution of the inner substruc- 
ture of the global jet close to the jet axis and the accretion disk. 
The size of this region is 60 x 20r^ and is comparable e.g. to the 
full grid in OP97. This part of the jet is not influenced from any 
outflow boundary condition or the jet evolution on the global 
scale. Note that the major part of the global jet is super fast 
magnetosonic, hence, there is no (physical) way to transport 
information from this part in upstream direction into the inner 
jet. 

4.3. Quasi-stationary nature of the (inner) jet 

As a major outcome of OP97 and FE00 it has been found that 
under a certain choice of boundary conditions the disk outflow 
may evolve into a stationary state after several hundreds of disk 
rotations. Such a stationary state solution provides an ideal test 
bed for investigation of the internal forces acting in the jet, 
since the well known equations of axisymmetric ideal MHD 
and its conservation laws may be used for interpretation of the 
jet flow. 

In our present simulations considering the MHD jet for- 
mation under the influence of magnetic diffusion, we also find 
such stationary states. We denote them as quasi- stationary 
since - due to the large computational domain - such a be- 
havior can be seen only in the inner region within reasonable 
computational time. The outer regions surrounding the (inner) 
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Fig. 6. Comparison of poloidal velocity vectors. Shown are the 
poloidal velocity vectors at the time of quasi-stationarity. Over- 
plot of the velocity field for the rj = simulation at t = 350 
with the ?] = 0.1 simulation at t = 400. The de-collimation 
visible in the poloidal magnetic field lines (see Fig. 5) is not 
present here. The thick line indicates the Alfven surface for 
r\ = 0.1, the thin line the fast magnetosonic surface. 



stationary jet flow will further evolve in time. For the area of 
the inner jet the stationary state is reached after approximately 
t = 350 disk rotations in the case of ideal MHD. This time 
scale can be quite different in the case of a non-vanishing mag- 
netic diffusion and also depends on the plasma-beta parameter 
fa. Note that in spite of considering magnetic diffusivity, we 
use the same parameter setup as OP97 and FE00. 

Figure 3 shows an over-plot of poloidal magnetic field 
lines resulting from our simulations at the time steps of t = 
250,300,350,400 disk rotations, respectively, in the case of 
relatively large magnetic diffusion rj = 0.1. It can clearly be 
seen how the poloidal magnetic field lines first diffuse outwards 
from their position at t = 250 (which is close to the non- 
diffusive field distribution, see also Fig. 5). After some hun- 
dred of rotations more, the field distribution approaches the 
quasi-stationary state (note the dotted lines almost matching 
the dashed ones). 

First, we note that the existence of such a quasi-stationary 
state might not necessarily be expected in the case of magnetic 
diffusion and instead one might think that the magnetic field 
will just decay forever. The reason, why a stationary state is 
possible in the case of a jet flow, is that with such a setup a loss 
of magnetic energy in the jet caused by diffusion can be replen- 
ished by the constant Poynting flux rising from the disk. This 
energy reservoir can be thought of as established constantly by 
the disk rotation and accretion and eventually by the gravita- 
tional potential of the central star. Note that also in the ideal 
MHD case the jet flow carries energy out of the computational 
box. Also this energy loss is balanced by the same way. The 
additional effect due to magnetic diffusivity is small compared 
to the total energy flow in the jet. 

In general, our simulations show that with increasing 77 the 
flow reaches the quasi-stationary state at a later time. For a 
smaller than moderate magnetic diffusivity (77 = 0, ...,0.5), we 



find an approximately linear relation between this time and the 
diffusivity (Fig. 4, right panel). This is also the range where the 
jet flow is a classical MHD jet like the Blandford-Payne solu- 
tion, starting as a sub-Alfvenic flow from the disk surface and 
crossing the Alfven surface at some height above the disk. For 
higher diffusivity the Alfven surface has been advanced into 
the disk for small radii and the character of the MHD flow is 
changed. We do not consider these cases in particular, however, 
we show the time of quasi-stationary state for completeness. 

For comparison, the local time step accounting for mag- 
netic diffusion and the Alfven time step, together with the 
global magnetic Reynolds number R m (defined by the jet size 
and a mean value for the jet Alfven speed, see Sect. 3.3) and 
the related a m are given in the Appendix in Tab. B.l for the 
time when the quasi-stationary state is reached. In Table B.l 
we have chosen a "typical" value for the Alfven speed within 
the grid of the inner jet. That this is feasible, is shown in Fig. 4 
(left) where we plot the variation of R m along the jet in the case 
77 = 0.1. As R m does not change along the jet by more than a 
factor of two (this is similar for other diffusivity), this value 
provides a good estimate for the global jet dynamical behavior. 

4.4. Jet velocity and collimation 

The most interesting, since directly observable, quantities of a 
jet are its velocity and degree of collimation. In Figure 6 we 
show the poloidal velocity vectors in the inner jet for a simu- 
lation without and with considering magnetic diffusion (77 = 
and rj = 0.1), taken at that time step, when the flow has reached 
the quasi-stationary state. 

The general point to mention is that the figures clearly show 
the self-collimating property of the MHD flow as the velocity 
vectors become more and more aligned with the jet axis as we 
go along the flow. However, there exists also a region of low 
collimation close to the disk where the velocity vectors point 
in radial direction (with ~ 45° half opening angle). 

For the low magnetic diffusivity (77 < 0.05) simulations 
we have observed an interesting feature. The apparent de- 
collimation of the poloidal magnetic field structure (Fig. 5), 
which is already present for a weak magnetic diffusivity (77 = 
0.01), is, however, not visible in the poloidal velocity (Fig. 6). 
In both cases, the flow evolution has reached a quasi-stationary 
state (t = 400). In contrast to ideal MHD, in the case of dif- 
fusive MHD, a miss-alignment between v p and B p is possible. 
Up to 77 = 0.5 the mismatch between the poloidal velocity and 
magnetic field vector is relatively weak for large z and about 
5° — 10° at intermediate heights above the disk. That means 
that we generally get a collimated stream along the axis, how- 
ever, the mass load distribution varies implying a variation of 
the mass flow rate through the r and z -boundaries with 77. If 
a larger magnetic diffusivity (7/ > 0.1) is applied, we observe 
also a de-collimation of the mass flow (see below). 

In general, both the poloidal and toroidal velocity increase 
with increasing diffusivity. At the same time, the magnetic field 
components decrease with the increasing diffusivity. This is 
shown in Fig. 7 where we plot the velocity and field compo- 
nents along jet at a distance of 15 r; from the jet axis. Thus, as 
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Fig. 7. Variation of the jet velocity and magnetic field with different magnetic diffusivity, 77 = (solid), rj = 0.1 (dashed), 
rj = 0.5 (dotted). Shown is the profile of the poloidal (left) and toroidal (right) components of the velocity (top) and the magnetic 
field (bottom) in z-direction at r = 15 at the time of quasi-stationarity. The velocity components increase, while the magnetic 
field strength decreases with increasing 77. Note the different scales for v p and v<f,. The boundary value for the toroidal velocity is 
the Keplerian value = 0.258 at r = 15. 



a conclusion, the diffusive jet becomes faster. We discuss that 
point below in the context of the Lorentz forces acting in the 
jet. 

However, two points should be mentioned concerning 
Fig. 7. The first is the decrease of toroidal velocity just above 
the disk surface. As the field line foot points rotate with 
Keplerian speed and the matter is corotating with the field lines, 
one would actually expect an increase of toroidal velocity, if 
the magnetic field is dominated by the poloidal component. 
In our case, we find at a certain radius along the disk surface 
that B^ = B p . Thus, the matter may substantially slide along 
the field. Just above the disk, the toroidal field strength first in- 
creases with height]], and the slide along the field line becomes 
larger. 

The second point is the fact that - unlike the poloidal veloc- 
ity - the magnetic field strength and the toroidal velocity do not 



2 This feature can be also observed in the simulations of OP97 (see 
their Fig. 4) 



seem to match the given boundary condition along the disk sur- 
face for the case of non-vanishing diffusivity. The reason is the 
jump in diffusivity between the disk (boundary condition) and 
the disk corona (active zones of the grid). The magnetic field 
lines are frozen-in the accretion disk, but when leaving the disk 
surface, they are immediately affected by magnetic diffusion. 
Thus, the magnetic field strength in the active zones of the grid 
(which are shown in Fig. 7) deviates from the boundary value 
in the case of a non vanishing diffusivity. In the stationary state 
solutions shown in Fig. 7 an equilibrium state has been reached 
between magnetic field diffusion and advection. We see that 
at this radius (r = 15) the field strength has increased for the 
region immediately above the disk. 



The observed de-collimation of the matter flow with in- 
creasing diffusivity is most evident if we plot the mass and mo- 
mentum fluxes across the boundaries of the inner jet region. We 
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define the fluxes across surfaces parallel to the accretion disk 
boundary by 



2nrpv z dr , Mv z = 



2-Krpvldr . (12) 



M 

lo Jo 
These are the mass flux and the kinetic z-momentum (i.e. mo- 
mentum in z-direction) flux along the jet axis. For the inner jet 
r mm = 20 and the integration is along z max = 60. The flux 
away from the jet axis (thus, in r-direction^|) is defined corre- 
spondingly by the integration along the r max -boundary from 
z = to z = z roa x. The corresponding flux into the jet (thus, 
prescribed by the disk boundary condition) has to be integrated 
along the 2 = axis. Signature of a good degree of collimation 
would be the fact that most of the mass and/or momentum flux 
is directed along the jet axis. 

Figure 8 shows how the mass and momentum fluxes for dif- 
ferent diffusivity evolve in time. We show the mass flux across 
the r and the z-boundary and the kinetic momentum flux in r 
and the z-direction integrated along both outflow boundaries. 

The large mass and momentum fluxes for the outflow dur- 
ing the first 100 — 200 rotations result from the fact that at these 
stages the initial hydrostatic corona is being pushed out of the 
grid of the inner jet. After the the bow-shock has left the inner 
grid, this initial coronal mass reservoir has been swept out and 
the remaining mass flow in the jet is given purely by the mass 
injection rate from the disk boundary condition. 

In the stationary state the mass inflow from the disk bound- 
ary into the jet must be equal to the mass loss across the bound- 
aries in r- and z-direction. That this is true in our simulations 
can be seen in Fig. 8 on the long time scale if we compare the 
solid line (inflow condition) to the sum of the dotted (radial out- 
flow) and the dashed lines (axial outflow). The analytical value 
for the mass rate as integrated from the given inflow bound- 
ary condition is M = 1.41 which is in good agreement with 
the numerical result. Note that for the momentum flux a simi- 
lar analytical integration gives a momentum flux from the disk 
into the jet of Mv z = 1.8 x 10~ 4 , which is much below the nu- 
merical value at the first active zone. As the momentum flux is 
not conserved as the matter becomes accelerated in the jet, this 
shows the tremendous gain of kinetic energy of the MHD flow. 
As a good estimate, the kinetic momentum flux in z-direction 
across the z = z max boundary is just the integrated mass flow 
rate M times the mean v z -velocity at this position. In the case 
of vanishing diffusivity we have Mv z ~ 1.5 x 0.6 = 0.9 which 
is similar to our numerically integrated momentum fluxf]. 

Now we compare the fluxes for simulations with differ- 
ent magnetic diffusivity. For the simulation run with rj = 0.5 
(Fig. 8, bottom panel) we have a mass inflow rate of 1.5 (in 
dimensionless units). The mass loss rate across the grid bound- 
aries is about 0.45 in z-direction and 1.05 in r-direction. 
Compared to the corresponding values in the case of 77 = 0, 
where about 70% of the mass flow leaves the box in z-direction, 
this clearly shows that the mass flow for rj = 0.5 is less colli- 
mated. This situation is even more evident for the simulations 



with higher diffusivity rj = 1.0, where, however, the mass in- 
jection from the disk boundary is partly super-Alfvenic (not 
shown). Thus, even if the velocity vectors have more or less 
the same direction for diffusivity up to rj = 0.5, the mass load 
along the stream lines is different due to the fact the matter, 
driven by centrifugal forces, may diffuse outwards across the 
magnetic field lines enhancing the mass flow rate in radial di- 
rection. 

At this point it might be interesting to discuss the results 
of recent diffusive MHD simulations of the jet formation out 
of the accretion disk (Kuwabara et al. 2000). These authors 



: Note that we define the momentum flux in z-direction across the 
rmax-boundary as Mv z = J max 2tyt pv z v r dr 

4 Almost no kinetic momentum flux in r-direction leaves the inner 
box across the z-outflow boundary 



find that the jet launching from the disk critically depends on 
the strength of the magnetic diffusivity. For small diffusivity, 
mass accretion in the disk and jet formation take place occa- 
sionally. For intermediate diffusivity the disk-jet system may 
reach a steady state. For high diffusivity the accretion rate and 
outflow rate decrease with diffusivity and may even vanish. In 
respect to our results, these results imply that the only way to 
launch a stationary MHD jet is indeed to allow for a reasonable 
amount of magnetic diffusion. Further, as we find less collima- 
tion for higher diffusivity or, equivalently, a weaker jet, such a 
state of stationary jet formation may become less important as 
the mass flow rates in the disk and the outflow decrease sub- 
stantially. It would therefore be of great importance to follow 
the simulations of Kuwabara et al. ( 2000 ) for longer time scales 
comparable to our runs. 

If we eventually define the degree of jet collimation by the 
mass flux across the jet boundaries, our simulations reveal the 
existence of a critical value of the magnetic diffusivity in this 
respect. In Fig. 9 the ratio of the mass flux leaving the grid in z- 
direction to that in r-direction is shown for different magnetic 
diffusivity for the time when the bow shock has left the inner 
box. This figure shows directly that for a high diffusivity mass 
flux ratio exceeds unity, indicating a weakly collimated mass 
flow. Figures 8 and 9 clearly show that for our model setup 
there exists a critical value of the magnetic diffusivity, rj CI . In 
the simulations with r\ < r/ cr the MHD flow evolves into a col- 
limated stream. In contrary, for 77 > rj CI the flow remains only 
weakly collimated. The actual value for the critical i] depends 
on the plasma beta /3j. In our standard setup chosen for Fig. 8 
and 9, we find r/ cr ~ 0.3, 

We note that the momentum flux gives somewhat differ- 
ent picture. The momentum flux in z-direction is always larger 
than that in r-direction. For our setup we obtain a ratio of 
about 5-8 when we compare the momentum fluxes in each 
direction. This demonstrates first the very high efficiency of 
rotating MHD flow in converting rotational kinetic energy into 
poloidal kinetic flux. In this respect, if we would define the de- 
gree of collimation by the momentum fluxes, our jets would 
perfectly collimated also for higher diffusivity. This leaves the 
question of how the degree of jet collimation is properly de- 
fined. Clearly, for diffusive MHD jets the field structure is not 
an accurate measure of propagation. What concerns the obser- 
vational appearance, the mass flow distribution (or actually the 
density distribution) would be the theoretical equivalent to the 
observed intensity (as long as no emission maps can be pro- 
vided by the simulations). 
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Fig. 9. Time evolution of the mass flow ratio between the ra- 
dial outflow boundary (mass flow in r-direction) and the axial 
outflow boundary (mass flow in z-direction) for different mag- 
netic diffusivity, 77 = 0,0.1,0.5 (solid, dashed and dotted line, 
respectively), in the inner part of the jet, (z X r) = (60 x 20)n. 
The final point of each line correspond to the end of the simu- 
lation when the (quasi-) stationary state has been reached. For 
higher diffusivity, the mass flux ratio in the quasi-stationary 
state increases indicating a decrease in degree of collimation. 



In summary, we propose that the mass flux gives the best 
measure of the degree of collimation. In our simulations we see 
a strong indication for the existence of a critical value of the 
magnetic diffusivity beyond which such a collimation cannot 
really be obtained. 

4.5. Lorentz forces in the jet 

Here, we deal with the question how the jet internal structure 
is modified by the effect of magnetic diffusivity as a result of 
our numerical simulations. Compared to the ideal MHD simu- 
lation (OP97) our results of (i) a de-collimation of the poloidal 
magnetic field structure for any value of magnetic diffusivity, 
(ii) a de-collimation of the hydrodynamic flow for strong dif- 
fusion and (iii) an increase of the jet velocity with increasing 
diffusivity are obtained purely by adding a physical magnetic 
diffusivity to the originally ideal MHD code. 

As the flow evolution in the MHD simulation results from 
a combination of various physical effects - magnetic and iner- 
tial forces, pressure and gravity - it is not straight forward to 
distinguish between these contributions. 

However, it seems to be clear that magnetic fields are the 
main driver for the flow acceleration and self-collimation and 
that, consequently, the addition of magnetic diffusion will mod- 
ify the MHD structure of the jet. Therefore, investigating the 



Lorentz forces in the quasi-stationary state may provide some 
insight into the physical mechanisms at work. 

At this point it might be instructive to recall in brief the 
basic mechanisms of MHD jet formation. Following the stan- 
dard model (e.g. Blandford & Payne 1982| , Ferreira 1997), a 
jet is launched as a sub-Alfvenic disk wind (by some unspeci- 
fied - but also unknown - process) and becomes accelerated 
by magneto-centrifugal forces in a strong poloidal magnetic 
field at first hand. As the flow approaches the Alfven surface, 
a toroidal magnetic field component is induced ("wound-up") 
due to the inertial back-reaction of the matter on the field. 
The toroidal field may lead to (de-) accelerating Lorentz forces 
-^Vh ~ J-L x B$ and (de-) collimating forces ~ Jn X B 
where, here, the perpendicular and parallel projection is made 
with respect to the poloidal magnetic field (which, only in the 
case of ideal stationary MHD is parallel to the poloidal veloc- 
ity). The toroidal Lorentz force F^^ ~ j p x B p affects the 
angular velocity of the matter, disturbing the centrifugal bal- 
ance and, thus, give rise also to a poloidal motion. 

Therefore, a change in the jet acceleration and collimation 
might be explained by the interplay of two mechanisms. First, 
the winding-up of the poloidal magnetic fields is less efficient 
since magnetic diffusion leads to a slip of matter across the 
field. As a consequence, the induced toroidal magnetic field is 
weaker leading to a less efficient acceleration by Lorentz forces 
but also to a de-collimation. This effect applies predominantly 
in the super-Alfvenic regime. Second, as a de-collimation of 
the poloidal magnetic structure also implies a smaller launch- 
ing angle for the sub-Alfvenic flow, the magneto-centrifugal 
acceleration mechanism may work more effective. As a conse- 
quence, the resulting fluid velocities in the jet should be larger, 
as indeed suggested by our simulations (Fig. 7). 

In Fig. 10 we show for different magnetic diffusivity the 
Lorentz force components along a field line (or, respectively, 
along the corresponding magnetic flux surface) leaving the nu- 
merical grid of the inner jet at (r = 20, z = 60). Note that due 
to the magnetic field de-collimation with 77, we compare differ- 
ent magnetic flux surfaces. These flux surfaces have their foot 
point between r = 5 and r = 8 along the disk surface and the 
Alfven point at about z = 25, 15, 5 (77 = 0, 0.1, 0.5). The fig- 
ure shows the Lorentz force components Fl.4> an d 
the corresponding acceleration of the fluid _Fl / P- 

The first point to mention is that the magnitude of the 
Lorentz force generally increases with increasing magnetic 
diffusivity. This is interesting insofar as the magnetic field 
strength decreases with increasing diffusion (Fig. 7). The 
Lorentz force has its maximum in that region before the Alfven 
point where the curvature of the poloidal field is largest. Thus, 
magnetic acceleration mainly works in this regime. 

This picture of an acceleration purely by magnetic forces 
is complementary to the above mentioned picture of an en- 
hanced magneto-centrifugal effect. For the parallel compo- 
nent this may directly lead to the observed increase in the 
poloidal velocity with increasing magnetic diffusivity (see 
Fig. 7). Additionally, the higher velocity also leads to stronger 
inertial forces and, for moderate heights above the disk, the dif- 
fusive plasma flow will tend to maintain its (radial) direction 
even if the field lines bend in direction of the jet axis. This will 
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Fig. 8. Time evolution of the mass flux and kinetic momentum flux for different magnetic diffusivity, r\ = 0, 0.1, 0.5 (top, middle 
and bottom figures, respectively), in the inner part of the jet, (z x r) = (60 x 20)r;. The final point of each line corresponds to 
the end of the simulation when the (quasi-) stationary state has been reached. Shown is the mass flux (left) across the different 
boundaries. The mass inflow across the first active grid cells along the (z = 0)-boundary (solid), across the outer (z = z max ) axial 
boundary (dashed), and across the outer (r = r max ) radial boundary (dotted). Also shown is the kinetic momentum flux across 
the boundaries (right). Poloidal kinetic momentum flux across the first active grid cells along the (z = 0)-boundary (solid). Note 
that this is already evolved from the value of the boundary condition. Momentum flux in z-direction integrated along the outflow 
boundaries (dashed). Momentum flux in r-direction integrated along the outflow boundaries (dotted). 
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Fig. 10. Lorentz forces in the jet for different magnetic diffusivity rj = 0, 0.1, 0.5 {solid, dashed, dotted lines). Left (Normalized) 
values of the force component perpendicular {top) and parallel {middle) to the field line and the toroidal component {bottom), 
Fl,±, F^n, Fl,4>, along a flux surface leaving the box of the inner jet close to {R = 20, Z = 60)-corner (see Fig. 3 or Fig. 5). 
For Fl.x the positive sign denotes the r-direction (de-collimating force). For F L u the positive sign denotes the z-direction 
(accelerating force). For Fl $ the positive sign denotes the (^-direction. Right Corresponding values of the magnetic acceleration 
(F Lt± /p),(F L>ll /p), (FlJp). 



re-distribute the mass flow distribution along the field line. The 
parallel component decreases rapidly with increasing z as it can 
be expected when the jet flow becomes more and more colli- 
mated. The same holds for the toroidal Lorentz force compo- 
nent. This component accelerates the plasma in toroidal direc- 



tion leading to an additional centrifugal effect which drives the 
matter in radial direction diffusing across the magnetic field. 
This is the reason for the increase of the mass flow rate along 
the outer stream lines with increasing diffusivity. 
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While the curves for the three Lorentz force components 
look quite similar at a first glimpse, we see that the correspond- 
ing components for the acceleration are somewhat different. 
The perpendicular (collimating) component^] of the accelera- 
tion remains on a rather high level throughout the (inner) box. 
That means that also in the asymptotic regime of the collimated 
jet these forces continue to collimate the jet flow. 

On the other hand, compared to the perpendicular compo- 
nent, the parallel and the (f> components of the corresponding 
acceleration have a steeper maximum and decrease to an only 
marginal strength beyond the Alfven point. This is what one 
would expect also from the standard MHD jet model. 

For 77 = and rj = 0.1, we see only a slight difference 
for the strength of the perpendicular components of force and 
acceleration at the large distances. Therefore, the degree of lo- 
cal flow collimation should be similar, as it is indeed visible 
in the poloidal velocity vectors, which are well aligned for the 
diffusivity considered (Fig. 6). However, we note the larger de- 
viation of the perpendicular components for 77 = 0.5, which 
mirrors the fact that in this case we are above the critical value 
rjci concerning the mass flow collimation (Fig. 8). 

In summary, our discussion of the Lorentz forces and its as- 
sociated acceleration gives a self-consistent picture of what we 
have observed in our numerical simulation. The perpendicular 
Lorentz force is essential for the collimation throughout the en- 
tire (inner) flow. The increase of the parallel Lorentz force for 
higher magnetic diffusivity gives rise to the higher velocities in 
the jet flow. The toroidal Lorentz force leads to an additional 
centrifugal effect enhancing the mass flow rate in the outer (yet 
un-collimated) parts of the jet flow. 

5. Summary 

In this paper we presented time-dependent simulations of the 
formation of axisymmetric protostellar MHD jets. In partic- 
ular, we were considering the effects of resistive MHD on 
the collimation and acceleration of the jet flow. Similar to re- 
cent si mulatio ns considering th e idea l MHD case (Ouyed & 
Pudritz 



1997 



Fendt & Elstner 2000), the accretion disk has 



been taken as a fixed boundary condition during the simulation, 
prescribing the mass flow rate and the magnetic flux distribu- 
tion. Our initial condition is a force-free magnetic field in a 
hydrostatic corona. Our simulations were performed on a grid 
of (z x r)=(280x 40) inner disk radii with 900 x 200 grid ele- 
ments or on a grid of 140 x 40 inner disk radii with 280 x 80 
grid elements. We find that in general the low resolution sim- 
ulations were sufficient to cover all physical effects observed 
in the higher resolution runs. In our discussion we mostly con- 
centrate on the structure of the inner jet which is the region of 
60 x 20 inner disk radii close to the star. We summarize our 
results as follows. 

(1) We have successfully implemented the physical magnetic 
diffusivity into the ZEUS-3D code. 



5 We note that the sign for -Fl,± is denned positive for the force 
vector pointing radially outwards. Thus, the increase of |-Fl.±| indi- 
cates an increase of the collimating Lorentz force on the matter. 



(2) We have discussed some analytical estimates about the 
strength of magnetic diffusivity in protostellar jets. We de- 
rive the distribution of magnetic diffusivity self-consistent 
to the turbulent Alfvenic pressure which is underlying our 
simulations and also present in the simulations of other au- 
thors. Our simulations, however, are performed with a con- 
stant diffusivity. Our main results do not depend on the ac- 
tual distribution of magnetic diffusivity. 

(3) In the global scale of our simulation, the jet bow shock ad- 
vances slower through the initial hydrostatic corona for the 
diffusive jets. The reason is the lower mass flux in the direc- 
tion along the jet axis in these jets. As expected, the internal 
structure of the jet is less disturbed in the case of diffusion. 
The Alfven surface comes closer to the disk surface. 

(3) For our model setup we find t hat, si milar to the case of idea l 
MHD jets (Ouyed & Pudritz |l997j Fendt & Elstner |2000| ), 
also resistive MHD jets can reach a quasi-stationary state. 
With increasing magnetic diffusivity, the quasi-stationary 
state of the jet is reached later. 

(4) With increasing diffusivity the jet velocity increases. The 
direction of the velocity vectors does, however, only change 
weakly. At the same time the poloidal magnetic field distri- 
bution becomes increasingly de-collimated. 

(5) As a proper measure of the degree of collimation we sug- 
gest the mass flux. If we compare the mass flow rates 
through the grid boundaries for different diffusivity we find 
strong indication for the existence of a critical value for 
the magnetic diffusivity r7 cr concerning the jet collimation. 
Beyond this value we still find an almost cylindrically col- 
limated stream along the jet axis, however, the bulk mass 
flow is in radial direction. For our setup, the critical (nor- 
malized) diffusivity is about rj CI w 0.3. 

(6) We discuss a self-consistent picture where these effects of 
jet de-collimation and acceleration are explained in the con- 
text of Lorentz forces. The perpendicular Lorentz force 
is essential for the collimation throughout the entire flow 
along the jet axis. The parallel Lorentz force increases for 
increasing magnetic diffusivity and gives rise to the higher 
velocities in the jet flow. The toroidal Lorentz force accel- 
erates the plasma in toroidal direction. This leads to addi- 
tional centrifugal forces re-distributing the mass flow rates 
across the magnetic flux surfaces towards the outer (yet un- 
collimated) parts of the flow. The latter two components 
play no role for larger distances along the flow. 



With our results we have shown that magnetic diffusivity 
plays indeed a role for the jet formation process. Turbulence as 
a natural (and necessary) property of accretion disks will natu- 
rally enter the disk wind and will be further advected into the 
jet. As we see in our simulations only a weak collimation for a 
high magnetic diffusivity, a hypothetical, and for sure exagger- 
ated, conclusion might therefore be that highly turbulent disks 
cannot drive a collimated jet mass flow. Such a claim may even- 
tually be tested by the astronomical observation and may also 
give some hint to answer the question why some disks have jets 
and some do not. 

Our present study should be understood as a first step in 
the right direction. Future work may improve the numerical 
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resolution and the grid size but may also consider e.g. an ad- 
ditional central stellar magnetosphere as a boundary condition. 
The most interesting (but also most difficult) prospect would be 
to include the evolution of the disk structure in the simulation. 

Appendix A: Numerical Tests 

As we have introduced the effect of physical magnetic diffusiv- 
ity into the ideal MHD ZEUS-3D code, careful tests were nec- 
essary to prove our implementation. In particular we checked 
the time scales introduced by magnetic diffusion and the behav- 
ior along the boundaries. The boundary conditions for an ax- 
isymmetric jet ("outflow", "inflow "and "reflecting") are quite 
different from what is e.g. needed in box simulations used for 
other scientific questions. 

We defined two numerical tests for our diffusive code. In 
both cases the code basically solves the diffusion equation. We 
obtained this limit by setting the initial density in the simulation 
to arbitrarily large values (here the normalized p ~ 10 9 ) effec- 
tively reducing any fluid motions in our simulations. The first 
test example is the analytical solution of the diffusion equation 
in Cartesian coordinates, the second example an axisymmetric 
torus of purely toroidal magnetic field in cylindrical coordi- 
nates. 

A. 1. Analytical solution to the diffusion equation 

In Cartesian coordinates (x, y, z) the solution of the one dimen- 
sional diffusion equation for infinite space is 



B z (y,t) = "^exp 



(y - Vof 
A-qt 



(A.l) 



7(47T<7) 



with the magnetic diffusivity 77 

As a test for our code, we choose as initial condition the 
magnetic field B z (x,y) = B z (y) for a certain time t — to from 
eq. (A.l). For the two-dimensional numerical grid we prescribe 
"free" (i.e. outflow) boundary conditions in ^-direction and a 
time-varying field for the boundaries in y-direction. 

Figure A. 1 show the result of our simulations for the time 
steps t = t At (At = 0.1,0.2,0.3,0.4) for a magnetic dif- 
fusivity 77 = 1.0 in comparison with the analytical results. As 
result, we obtain a perfect agreement between the numerical 
simulation and the analytical solution. 

A. 2. Toroidal field torus 

Here, our aim is to check how our code treats magnetic diffu- 
sion in cylindrical coordinates, along the outflow boundary in 
r-direction and along the symmetry axis. As initial condition, 
we define a torus of toroidal magnetic field 



B v {r,z,t Q = 0.1) = ^-exp | 



(r - r ) 2 + (2 - z ) : 



Figure A. 2 show the result of our simulations for the time steps 
t = t At (At = 0.1,0.2,0.3,0.4) for a magnetic diffusiv- 
ity 77 = 1.0 in comparison with the analytical results. The 
simulation shows how the peak of the field distribution moves 
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Fig. A.l. Numerical test of magnetic diffusion. Grid size 50 x 
50 elements for a normalized physical grid 2.0 x 2.0. Left 
Isocontours of the magnetic field strength B z (x,y) (normal- 
ized units) for different time steps t = to + At, At = 0.0 (solid 
line), At = 0.1 (dashed line). Right Normalized intensity pro- 
file of the magnetic field strength across the two-dimensional 
box along x = 1 for different time steps t = to + At with 
At = 0.0, 0.1, 0.2, 0.3 (top to bottom curve). Comparison be- 
tween the analytical solution solid lines and the numerical sim- 
ulation (dashed lines). 
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Fig. A.2. Numerical test of magnetic diffusion in cylindrical 
coordinates. Grid size 100 x 50 elements for a normalized phys- 
ical grid 10.0 x 5.0. Left Isocontours of the toroidal magnetic 
field strength for different time steps t = to + At, At = 0.0 
(solid line), At = 0.1 (dashed line). Right Normalized in- 
tensity profile of the magnetic field strength across the two- 
dimensional box along z = 5 for different time steps t = 
t + At with At = 0.0, 0.1, 0.2, 0.3 (top to bottom curve). 

slightly inwards from its initial central position as the field dif- 
fuses. Note that as the field diffuses outwards the volume over 
which the toroidal field is distributed increases. Therefore the 
decrease in field strength in outward direction. Along the sym- 
metry axis the field strength remains zero, whereas the field 
strength along the outflow boundaries increases. No boundary 
condition is prescribed here. 

Although there is no analytical solution to compare with, 
this simulation gives again convincing evidence that we prop- 
erly incorporated the magnetic diffusion in the ZEUS-3D code. 

Appendix B: Time scales in the simulation 

Table 1 shows the time scales and typical number values ap- 
plied in our jet simulations for different magnetic diffusivity 
77. As an example we refer to the simulation runs with rela- 
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Table B.l. Typical numbers for our simulations for different 
magnetic diffusivity 77. Diffusive time step t v , Alfven time step 
Ta, the global magnetic Reynolds number R m and the mag- 
netic turbulence parameter a m . 



V 


r„ 


TA 


Rm 


a m 





00 


0.258 


OO 





0.1 


0.625 


0.293 


40 


0.025 


0.15 


0.417 


0.084 


30 


0.033 


0.25 


0.25 


0.136 


20 


0.05 


0.5 


0.125 


0.080 


10 


0.1 


1 


0.063 


0.183 


8 


0.125 


2.5 


0.025 


1.394 


4 


0.25 



tively low numerical resolution. The diffusive time step (i.e. 
the local time scale for magnetic diffusion) is T n = (As/77 2 ) 
with the size of the grid element Ax. For comparison we show 
the Alfven time step ta = (Ax/va)- The other two numbers 
consider the global evolution of the jet flow. For the magnetic 
Reynolds number R m = (v\r max /rf) we have considered the 
global size of the jet and similar for the magnetic turbulence 
parameter a m . 
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